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An  experimental  investigation  is  carried  out  to  characterize  the  performance  of  thermoelectric  modules 
used  for  electric  power  generation  over  a  range  of  different  resistance  loads.  The  performance  of  a  Peltier 
cell  used  as  a  thermoelectric  generator  is  evaluated  in  terms  of  power  output  and  conversion  efficiency. 
The  results  show  that  a  thermoelectric  module  is  a  promising  device  for  waste  heat  recovery. 
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1.  Introduction 

Thermoelectric  conversion  may  be  defined  as  the  result  of  a 
process  by  which  heat  is  converted  into  electricity  through  the 
use  of  a  heat-to-electricity  conversion  device.  Thermoelectric  con¬ 
version.  besides  occurring  in  simple  thermocouples,  forms  the  ba¬ 
sis  of  reversible  thermoelectric  modules  (TEMs),  which  can  work 
either  as  coolers  when  supplied  with  electricity,  or  as  generators 
when  supplied  with  heat. 

A  TEM  consists  of  alternate  ingot-shaped  n-  and  p-type  semi¬ 
conductor  thermoelements,  which  are  connected  electrically  in 
series  with  metal  connecting  strips,  sandwiched  between  two  elec¬ 
trically  insulating  but  thermally  conducting  ceramic  plates  [1], 

The  utilization  (generator  or  cooler)  determines  the  tempera¬ 
ture  under  which  the  thermoelements  operate,  and  hence  the 
choice  of  materials  for  the  thermoelements.  For  low  temperature 
applications  (commonly  coolers)  BiTe  alloys  are  used;  in  contrast 
high  temperature  materials  such  as  PbTe  or  SiGe  alloys,  are  used 
in  thermoelectric  generators  [2],  However,  the  BiTe  alloys  used  in 
conventional  coolers  have  recently  become  relatively  cheap,  sug¬ 
gesting  their  use  also  in  thermoelectric  generators  [3]. 

Recent  review  papers  (Riffat  and  Ma  [4],  Bell  [5],  Kajikwa  [6], 
and  for  the  automotive  sector  Vazquez  et  al.  [7])  summarize  the 
potential  applications  in  power  generation  of  thermoelectric 
devices.  These  works  make  it  evident  that  recent  advances  in 
materials  and  material  processing  have  led  to  higher  theoretical 
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thermoelectric  conversion  efficiencies.  However,  when  using  TEMs 
in  the  power  generation  mode  the  final  thermal  efficiency  is  lim¬ 
ited  also  by  other  factors:  the  reduction  of  the  heat  losses,  the  de¬ 
sign  optimization  to  reduce  the  parasitic  losses  and  the  appropriate 
choice  of  materials. 

In  the  literature  only  a  limited  number  of  papers  are  available 
dealing  with  the  thermal  experimental  performance  of  a  thermo¬ 
electric  generator.  Some  authors  investigate  directly  the  perfor¬ 
mance  of  TEMs.  Leavitt  et  al.  [8]  discuss  the  critical  requirements 
for  accurate  testing  of  a  TEM  and  Takazawa  et  al.  [9]  report  the  de¬ 
sign  and  development  of  a  testing  system  for  TEMs  operating  at 
high  temperature.  Rauscher  et  al.  [10]  present  an  apparatus  for 
measuring  the  conversion  efficiency  and  other  key  properties  of 
TEMs;  experimental  data  for  a  high  efficiency  Bi2Te3  module  are 
discussed  and  compared  with  the  results  obtained  from  a  theoret¬ 
ical  model.  From  their  experimental  data  Sandoz-Rosado  and  Ste¬ 
vens  [11]  analyze  a  TEM  under  a  wide  range  of  temperature  and 
loading  conditions.  Their  attention  focused  on  the  measurement 
of  the  heat  flow  into  the  module  and  the  parasitic  losses  inside 
the  TEM.  Gou  et  al.  [12]  study  experimentally  a  low- temperature 
waste  heat  thermoelectric  generator.  The  data  are  used  for  the  val¬ 
idation  of  the  model  used  to  guide  the  analysis  and  the  optimiza¬ 
tion  of  the  generator. 

Other  authors  carry  out  direct  investigations  on  specific  appli¬ 
cations  of  TEMs.  Crane  and  Jackson  [13]  study  the  behaviour  of  a 
thermoelectric  generator  heated  by  a  hot  liquid  flow  and  cooled 
by  an  air  stream.  Saqr  et  al.  [14]  report  the  application  of  an  ex¬ 
haust-based  thermoelectric  generator  in  the  thermal  design  of  an 
automobile.  For  three  different  examples  they  obtained  data  for 
the  overall  efficiency,  the  temperature  difference  across  the  TEMs 
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Nomenclature 

A  cross  area  of  the  thermoelement  (mm2) 

& 

thermal  grease 

I  electric  current  (A) 

gen 

thermoelectric  generator 

L  length  (mm) 

H 

hot  junction 

K  thermal  conductance  (W  I<  ') 

In 

input 

m  number  of  modules  in  the  generator 

ins 

insulation 

N  number  of  thermocouples  per  module 

J 

Joule  heater 

P  electric  power  output  on  the  load  resistance  (W) 

L 

resistance  load 

Q,  heat  power  (W) 

n 

n-type  semiconductor 

R  electric  resistance  (£1) 

0  V 

overall 

T  temperature  (K) 

P 

p-type  semiconductor 

V  electric  voltage  (V) 

TEM 

thermo  electric  module 

Z  figure  of  merit  (ot2/p2)  (K  ') 

OC 

open  circuit 

s 

stainless  steel 

Subscripts 

A  air 

Greek  symbols 

b  bolts 

a 

Seebeck  coefficient  (VK  ') 

CC  closed  circuit 

A 

difference 

C  cold  junction 

X 

thermal  conductivity  (W m  1  K  ') 

cer  ceramic  insulator 

P 

electrical  resistivity  (£2  m) 

env  environment 

e 

conversion  efficiency 

and  the  conversion  efficiency.  Niu  et  al.  [15]  carry  out  an  experi¬ 
mental  study  of  a  thermoelectric  conversion  unit,  consisting  of 
commercially  available  TEMs  incorporated  in  a  parallel  plate  heat 
exchanger. 

In  the  process  of  reduction  of  the  direct  experimental  data,  it  is 
necessary  to  use  a  mathematical  model  of  greater  or  less  simplifi¬ 
cation.  Such  a  model  involves  the  electrical  and  thermal  properties 
of  the  TEM,  and  often  these  are  unknown.  Interesting  procedures 
have  been  recently  developed  to  obtain  these  properties  from 
manufacturers’  data  (Luo  [16]  and  Palacios  et  al.  [17]).  These  are 
preferable  to  the  situation  where  the  experimental  data  them¬ 
selves  are  used  to  determine  the  parameters  [13]. 

The  present  paper  reports  an  experimental  investigation  of  the 
performance  of  a  power  generation  device  in  which  we  have  used 
multiple  Peltier  modules  in  the  Seebeck  mode.  The  thermoelectric 
generator  is  analyzed  on  the  basis  of  the  experimental  data  for  the 
‘open’  and  ‘closed’  circuit  voltage,  electric  power  output  and  con¬ 
version  efficiency  as  a  function  of  the  temperature.  In  the  labora¬ 
tory  a  Joule  resistor  is  used  as  the  thermal  source  instead  of 
waste  heat.  The  experimental  data  give  significant  information 
about  the  behaviour  of  the  thermoelectric  generator.  In  particular 
the  effects  of  practical  but  unavoidable  details  are  evidenced,  such 
as  the  clamping  of  bolts  and  the  fitting  of  thermal  insulation  be¬ 
tween  TEMs. 


2.  Experimental  apparatus 

A  schematic  diagram  of  the  test  section  is  shown  in  Fig.  1.  The 
test  section  consists  of  a  thermoelectric  generator  made  of  nine 
TEMs  (40  x  40  x  3.8  mm)  connected  electrically  in  series.  They 
are  sandwiched  between  a  square  plate  of  aluminium  (200  x 
200  x  4  mm)  and  a  conventional  heat  sink  of  the  same  area. 
Uniform  pressure  is  ensured  by  using  22  stainless  steel  screws 
(2.5  mm  nominal  diameter,  Xs  =16Wm_1r')  to  clamp  the 
assembly  together.  Uniform  contact  between  the  Peltier  cells  and 
the  aluminium  plates  is  enhanced  by  a  thin  layer  of  thermal  grease 
(SG500,  /g=  0.77  W  m-1  K  ’).  Thermal  insulation  comprised  of  cel¬ 
lular  rubber  is  placed  between  the  modules. 

The  overall  dimensions  of  the  aluminium  heat  sink  are  200  mm 
by  39.2  mm  by  200  mm  for  the  width,  height  and  length  respec¬ 


tively.  The  finned  surface  consists  of  20  fins  with  thickness 
3.2  mm  and  height  30.8  mm.  The  gap  between  the  fins  is  7.2  mm. 
To  facilitate  heat  release  by  the  thermoelectric  generator,  the  fins 
of  the  bottom  plate  are  immersed  in  a  bath,  where  water  flows 
at  room  temperature  (Photo  1). 

The  hydraulic  circuit  is  schematically  shown  in  Fig.  2.  It  is  of  the 
‘open’  type,  that  is  with  constant  head  feeding  and  receiving  tanks 
open  to  the  atmosphere.  The  former  feeds  the  heat  transfer  section. 
After  passing  this  section,  the  fluid  is  collected  into  the  receiving 
tank,  from  where  it  is  pumped  up  to  the  feeding  tank.  The  operat¬ 
ing  fluid  is  cooled  through  the  pipe  walls  by  heat  transfer  with  the 
surrounding  air. 

The  measured  flow  rate  of  the  water  is  7.2  x  10  s  m3  s  \  In  the 
gap  between  the  fins  the  Reynolds  number  is  272.  The  overall 
thermal  resistance  of  the  heat  sink,  based  on  the  experimental 
data,  is  0.0152  ±  0.0015  K/W.  An  estimation  of  the  overall  thermal 


Fig.  1.  Schematic  diagram  of  the  test  section. 
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Photo  1.  Thermoelectric  generator  in  a  small  experimental  bath. 


Table  1 

Geometrical  parameters  of  the  generator. 


Data  Symbol 

Cellular  rubber  thickness  lins 

Hot  side  thermal  grease  thickness  LgH 

Cold  side  thermal  grease  thickness  IgC 

Diameter  of  the  bolts  Db 

Average  thickness  of  the  heater  Theater 

Area  of  the  heater  Aheater 


The  electric  resistance  of  the  Joule  heater  has  been  measured  for 
different  values  of  the  temperature.  A  variation  law  of  this  electri¬ 
cal  resistance  is  obtained  by  means  of  a  linear  interpolation.  The 
knowledge  of  the  dependency  of  the  electric  resistance  on  the  tem¬ 
perature  allows  the  monitoring  of  the  electric  power  supplied  by 
means  of  simple  measurements  of  voltage. 

3.  System  modelling 


receiving  tank 


Fig.  2.  Schematic  overview  of  the  heat  transfer  system. 


resistance  of  the  heat  sink,  based  on  the  results  of  Morini  [18], 
gives  0.0148  K/W,  in  very  good  accordance. 

A  heat  flow  is  supplied  to  the  top  surface  of  the  thermoelectric 
generator.  It  is  obtained  by  dissipating  via  the  Joule  effect  an  as¬ 
signed  power  in  a  nickel-chrome  resistor.  This  resistor  of  square 
shape  and  of  the  same  dimensions  as  the  outside  top  surface  of 
the  thermoelectric  generator,  is  placed  under  a  layer  of  insulating 
material  (mineral  wool  near  the  heater,  glass  wool  near  outside). 
The  total  thickness  of  the  layers  of  insulating  material  is  200  mm 
and  the  thermal  conductivity  is  0.03  W  m  1  K  '.  The  insulating 
layer  is  confined  sideways  by  a  40  mm  thick  layer  of  expanded 
polyurethane. 

In  Table  1  are  reported  the  remaining  geometrical  parameters  of 
the  generator. 

The  temperature  on  the  surface  of  the  modules  is  measured 
using  copper-constantan  thermocouples.  Three  thermocouples 
are  placed  on  the  aluminium  finned  plate  in  the  space  between 
the  modules.  Nine  thermocouples  are  fixed  on  the  top  plate  over 
the  modules.  Two  armored  thermocouples  of  type  K  are  placed  in¬ 
side  the  duct  feeding  the  bath  and  in  the  air,  respectively. 

The  fourteen  thermocouples  are  connected  to  an  acquisition 
data  system  Agilent  34970A.  The  same  instrument  is  used  for  the 
measurement  of  the  circuit  voltage  on  the  resistance  load.  When 
more  thermocouples  are  placed  on  a  surface,  the  average  value  of 
the  data  is  considered. 


In  order  to  predict  the  performance  of  the  thermoelectric  gener¬ 
ator  as  a  function  of  the  experimental  thermal  conditions,  a 
lumped  parameter  model  has  been  formulated.  The  electrical 
scheme  of  the  model  is  shown  in  Fig.  3.  In  Fig.  4  a  section  of  a  ther¬ 
moelectric  module  is  shown. 

The  theoretical  model  gives  the  values  of  the  main  parameters, 
expressed  in  terms  of  the  temperature  difference  established  by 
the  heat  flow. 

Each  pair  of  thermoelements  (Fig.  4)  produces  an  electromotive 
force  by  the  Seebeck  effect  given  by: 

Voc  =(ocp-ocn)(T'H-T'c)  (1) 


Heat 

Source 


open  /  closed 


Fig.  3.  Electrical  scheme. 


Fig.  4.  Thermoelectric  i 
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The  top  and  bottom  temperatures,  T'H  and  T'c  respectively,  are  at 
the  hot  and  cold  junctions  of  the  thermoelements;  ap  and  y.„  are 
the  Seebeck  coefficients  of  the  p-  and  n-type  semiconductors, 
respectively. 

When  a  resistance  load  is  inserted,  the  power  production  starts. 
The  ‘closed’  circuit  voltage,  VCc,  at  the  ends  of  the  resistance  load, 
Rl,  is  given  by: 

Vcc  =  Rr£+RimN{a”  ~  an)(J'H  ~  r'c)  (2) 


where  the  electromotive  force  is  extended  to  the  number  of  ther¬ 
mocouples  per  module,  N,  and  the  number  of  modules,  m.  RTEM  is 
the  internal  resistance  of  the  thermoelectric  generator,  calculated 
as  the  sum  of  the  electric  resistance  of  each  thermoelement  pair,  ex¬ 
tended  as  above  to  the  numbers  N  and  m,  where  the  electric  resis¬ 
tance  per  thermoelement  pair  is  given  by  the  sum  of  the  electric 
resistance  of  two  p-  and  n-type  thermoelements: 


R'  =  Pnj 


(3) 


In  Eq.  (3)AP,  Lp,  pp  and  An,  Ln,  pn  indicate  the  cross  sectional  area,  the 
length  and  the  electrical  resistivity  of  the  p-  and  n-type  thermoele¬ 
ments,  respectively. 

The  power  output  dissipated  by  the  resistance  load  is  given  by: 


P  = 


Rl 

(Rtem  +  Rl)2 


[mN(ap  -  a „)(T'H  -  TE)]2 


(4) 


The  heat  input,  Qin,  for  the  pair  of  thermoelements  for  steady 
operations,  is  composed  of  three  terms  due  to  heat  conduction, 
and  the  Peltier  and  Joule  effects.  This  heat  flow  is  given  by: 

Qin  =  K'AT'  +  (Op  -  <xn)T'HI  -  ^R72  (5) 


where  K'  is  the  thermal  conductance  of  a  thermoelement  pair,  given 
by: 


In  Eq.  (6)  Ap,  lp  and  Am  indicate  the  cross  sectional  area  and  the 
thermal  conductivity  of  the  p-  and  n-type  thermoelements, 
respectively. 

The  heat  input  into  the  thermoelectric  modules  is  given  by: 
Qtem  =  mNQin  (7) 

where  N  and  m  are  as  defined  above. 

The  electric  current,  /,  in  Eq.  (5)  is  that  circulating  in  the  gener¬ 
ator  when  the  load  is  inserted,  given  by: 

mN(ap-an)(rH-rc) 

Rtem  +  Rl 

The  thermoelements  are  sandwiched  between  two  ceramic 
plates  (Fig.  4)  in  thermal  contact  with  the  assembly  surfaces 
through  thin  layers  of  thermal  grease.  The  temperature  is  mea¬ 
sured  on  the  assembly  surfaces,  firstly  in  the  immediate  proximity 
of  the  TEM  on  the  bottom,  and  also  on  the  top  plate  over  the  mod¬ 
ules  on  the  top.  Then,  to  obtain  the  values  of  temperature  inside 
the  ceramic  plates  at  the  hot  and  cold  junctions,  the  theoretical 
model  proposed  by  Chen  et  al.  [19]  is  used.  This  mathematical 
model  is  based  on  a  thermal  balance  for  the  TEM  and  its  connecting 
structure.  For  the  closed  circuit  the  temperature  difference  over 
the  junctions  ( T'H  -  T'c)  is  a  function  of  that  outside  the  ceramic 
plates  (Th  -  Tc)  and  other  parameters: 

(T'h  ~  T'c )  =  (Th  -  Tc)f(RL ,  Kh ,  Kc)  (9) 

where  RL  is  the  resistance  load,  Kh,  Kc  are  the  upstream  and  down¬ 
stream  thermal  conductance  of  the  hot  and  cold  junction  of  the  pair 
of  thermoelements,  respectively. 


When  the  circulating  current,  /,  tends  to  zero,  the  relationship 
between  AT  and  AT  becomes  (Rowe  [20]): 


r  -  ,  tH  lc  ' 


(10) 


In  this  case  for  each  thermoelement  pair  the  voltage  of  the  model 
tends  to  the  ‘open’  circuit  voltage  and  is  obtained  as: 

V'oc  =  (ap-an)  {Th~Tc)  (11) 

(1+£  +  fe) 

For  the  generator  as  a  whole,  and  using  the  data  from  all  the  pairs  of 
thermoelements,  the  ‘open’  circuit  voltage  Voc  is  given  by: 

Voc  =  mNV'0C  (12) 


4.  Experimental  tests 


Each  measurement  was  carried  out  with  the  load  in  (‘closed’  cir¬ 
cuit)  and  then  repeated  with  the  load  out  (‘open’  circuit).  To  start 
any  test  series  the  resistance  load  is  connected.  The  test  continues 
with  the  insertion  of  the  hydraulic  pump  to  circulate  the  water  in 
the  cooling  system.  Successively,  the  Joule  heater  is  supplied  to 
give  the  assigned  power  input.  The  monitoring  of  the  test  deter¬ 
mines  when  steady  state  conditions  have  been  attained,  which  is 
when  the  data  acquisition  is  started.  All  the  parameters  are  then 
sampled  every  150  s,  over  a  period  of  40  min. 

When  the  data  acquisition  is  finished  for  the  ‘closed’  circuit,  the 
resistance  load  is  disconnected  and  the  procedure  is  repeated  for 
the  ‘open’  circuit. 

The  useful  power  dissipated  by  the  resistance  load,  RL,  is  calcu¬ 
lated  on  the  basis  of  the  electric  measurements  using  the  following 
equation: 


where  VCc  is  the  ‘closed’  circuit  voltage  over  the  resistance  load. 

In  order  to  obtain  the  efficiency  of  the  generator  this  useful 
power  has  to  be  compared  with  the  heat  flow  input.  A  lumped 
parameter  distribution  of  the  heat  flow  can  be  obtained  with  a 
thermal  resistance  scheme  (Fig.  5).  The  overall  heat  power  input 
into  the  generator,  Qgen,  is  the  difference  between  the  electric 
power  supplied  to  the  Joule  heater,  Qj  and  the  heat  leakage  through 
the  insulation  materials  placed  over  the  Joule  heater,  Qenv: 

Qgen  =Qj-  Qenv  (14) 


Qgen  is  composed  of  two  terms:  the  heat  power  input  into  the 
TEMs,  Qtem.  and  the  heat  lost  by  conduction  through  the  assembly 
structures.  Consequently,  the  heat  power  into  the  TEMs  is  given 
by: 

Qtem  =  Qgen  “  (Qb  +  Qins)  (15) 


where  Qh  is  the  heat  leakage  through  the  bolts  clamping  the  assem¬ 
bly  together  and  Qins  is  the  heat  flow  through  the  insulating  mate¬ 
rial  (cellular  rubber)  placed  between  the  cells. 

Owing  to  the  significant  temperature  difference  between  the 
plates,  in  the  calculation  of  the  heat  losses  through  the  insulating 
material  a  linear  dependence  of  the  thermal  conductivity  on  tem¬ 
perature  is  adopted  for  the  cellular  rubber. 

The  conversion  efficiency  of  the  TEMs  is  given  by  the  ratio  be¬ 
tween  the  measured  power  output,  Eq.  (13),  and  the  net  power  into 
the  modules,  Eq.  (15): 


(16) 
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Heat 

Source 


water 


Fig.  5.  Thermal  scheme. 


The  overall  conversion  efficiency  of  the  generator  is  given  by 
the  ratio  between  the  measured  power  output,  Eq.  (13),  and  the 
overall  power  into  the  generator,  Eq.  (14): 


5.  Results  and  discussion 

Data  for  the  experimental  performance  of  the  generator  for  dif¬ 
ferent  operating  conditions  are  presented. 


In  this  experiment  very  cheap  commercially  available  Peltier 
cells  are  used  but  the  manufacturer’s  data  sheets  were  unavailable. 
Thus  the  procedures  proposed  in  [16,17]  were  inapplicable.  Also  the 
way  of  using  the  experimental  data  themselves  to  determine  the 
parameters,  as  adopted  in  [13],  was  considered  inappropriate.  Since 
many  parameters  are  mandatory  when  using  the  lumped  parameter 
model  for  predictions,  a  cell  was  destroyed  to  characterize  its 
components.  These  were  measured  using  an  optical  microscope. 
Referring  to  Fig.  4  for  the  meaning  of  the  symbols,  L  =  1.6  mm,  Lc 
=  0.5  mm,  Lcer  =  0.6  mm;  the  p  and  n  legs  show  the  same  square 
section  with  side  1.5  mm  and  transverse  area  A  =  2.25  mm2.  Every 
module  consists  of  127  pairs  of  thermoelements. 

With  X-ray  spectroscopy  the  chemical  composition  of  the 
semiconductors  was  identified:  Bio.5Sb1.5Te3  for  the  p-type  semi¬ 
conductor  and  Bi2Te3  for  the  n-type  semiconductor. 

For  the  remaining  parameters  which  are  needed,  values  are 
used  as  found  in  the  literature.  In  Tables  2  and  3  are  reported  the 
electrical  and  thermal  properties  of  the  two  semiconductors 
[21,22],  and  the  remaining  generator  parameters,  respectively. 

The  tests  were  intended  to  have  the  same  constant  temperature 
on  the  bottom  surface.  However,  due  to  the  cooling  system,  using 
the  ambient  air  as  the  sink,  this  constant  value  was  not  attained 
(32.7  ±  3.9  °C,  Tables  4  and  5).  The  variation  on  the  bottom  surface, 
even  if  limited  due  to  the  efficient  heat  transfer  between  air  and 
the  circulating  water,  is  mainly  due  to  variation  in  the  ambient 
temperature. 

When  comparing  the  experimental  data  with  the  results  of  the 
model,  instead  of  a  point  to  point  comparison  between  data  and 
predictions,  a  graphical  comparison  for  groups  of  data  (typically 
at  the  same  power  input  or  at  the  same  resistance  load)  is  adopted. 
To  make  this  comparison  clear,  in  the  predictions  a  single  reference 
temperature  was  chosen  for  the  cold  surface.  This  temperature  is 
the  average  of  those  measured  (32.7  °C,  Tables  4  and  5). 

Unlike  typical  data  reductions  which  report  voltage  and  power 
produced  as  a  function  of  the  circulating  current  ([6,11]),  in  this 
work  particular  emphasis  is  laid  on  the  direct  dependency  of  the 
parameters  on  the  thermal  quantities  of  the  problem. 


Table  2 

Semiconductor  properties  referred  to  T  =  323  K. 

Semiconductors  Compound  a(pVK_1)  p(pfim)  2(Wm~'K~’)  Z(K~')  References 

p-type  Bio.sSb,  5Te3  130  5.26  1.69  1.9  x  10  3  [21] 

n-type  Bi2Te3  -213  26.10  0.823  2.11  x  103  [22] 


Table  3 

Thermophysical  properties  of  the  materials. 


Data  Symbol 


Thermal  conductivity  of  the  ceramic  plates  [23]  2cer 

Thermal  conductivity  of  the  cellular  rubber  [25]  Ai„s 

Thermal  conductivity  of  the  aluminium  [26]  2ai 

Thermal  conductivity  of  the  bolts  [26] 

Thermal  conductivity  of  the  heater  [26]  tNiCr 


Value  (W  m-1  IT1) 


36 

0.036  +  10~4T  (°C) 

194 

16 

16 


Table  4 

Top  and  bottom  experimental  temperature  for  the  open  circuit  condition. 
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Table  5 

Top  and  bottom  experimental  temperature  for  the  closed  circuit  condition. 


Qj  Ki  =  1 


11.14  Q 


R3  =  39.20  Cl 


R4  =  86.22  Cl 


Th  Tc  Ta  Tj  Th  Tc  Ta  Tj  Th  Tc  Ta 


Th  Tc  Ta  Tj 


108.4  55.0  30.1 

130.5  60.0  30.0 

167.1  71.7  33.5 

210.4  81.5  34.1 

240.7  87.1  33.5 


21.6  162.0  56.3 

21.2  186.6  61.5 

22.9  236.8  72.3 

21.8  288.4  83.8 

19.4  320.1  89.1 


30.2  22.7  163.3 

30.6  21.5  188.3 

33.2  21.7  236.9 

34.9  22.7  290.7 

34.2  19.9  321.4 


55.1  28.9  21.1 

61.4  29.9  20.8 

71.1  31.5  20.3 

83.6  34.5  22.1 

91.6  35.1  21.3 


158.2  56.6  29.4 

185.0  62.1  29.5 

231.5  73.5  31.8 

284.8  85.2  34.2 

322.6  94.4  36.0 


21.4  159 

19.7  184 

21.1  234 

21.4  283 

21.6  320 


0  250  500  750  1000  1250  1500  1750  2000  2250 

time  (s) 


Fig.  6.  Temperature  distribution  in  a  test  (closed/open  circuit). 


The  measurements  are  carried  out  for  five  values  of  the  electric 
power  supplied  to  the  heater,  Q/I  for  each  value  of  Qj  four  different 
resistance  loads  are  changed. 

In  Fig.  6  typical  temperature  distributions  are  shown  for  water, 
air,  bottom  and  top  surfaces  of  the  modules  and  top  surface  of  the 
insulating  layer,  for  the  cases  of  ‘closed’  and  ‘open’  circuit,  respec¬ 
tively.  The  two  test  cases  are  held  at  the  same  conditions.  The  tem¬ 
perature  difference  between  the  bottom  and  top  surfaces  of  the 
modules  is  greater  in  the  case  of  the  ‘open’  circuit  than  for  the 
‘closed’  circuit.  In  both  cases  the  bottom  temperature  is  about 
31  °C,  while  the  top  temperature  is  71  °C  and  75  °C,  for  the  ‘closed’ 
and  ‘open’  circuit,  respectively. 

In  Tables  4  and  5  the  temperatures  are  reported  for  the  TEMs 
surfaces,  heater  top  surface  and  air  for  the  whole  experimental 


set  of  trials.  On  the  bottom  surface  the  temperature  is  quite  con¬ 
stant  (32.7  ±  3.9  °C),  while  on  the  top  the  temperature  varies  from 
55.0  °C  to  98.3  °C.  The  spread  of  the  top  surface  temperature  de¬ 
pends  on  the  heat  power  input  supplied  (Tables  6  and  7). 

In  Tables  4  and  5,  for  a  given  electric  power  supplied,  Qj,  and  dif¬ 
ferent  resistance  loads,  the  top  surface  temperature  of  the  Joule 
heater,  7},  does  not  change  significantly.  Since  the  heat  leakage  to 
the  environment,  Qsm,  depends  on  7}  and  on  the  ambient  temper¬ 
ature,  the  heat  power  input  into  the  thermoelectric  generator,  Qgen 
in  Eq.  (14),  seems  to  be  independent  of  the  resistance  load. 

In  Tables  6  and  7  are  summarized  the  values  of  Qgen  and  QrEM 
obtained  from  the  measured  values  of  Qj  with  the  procedure  de¬ 
scribed  in  Section  4.  Since  both  the  heat  power  input  into  the  ther¬ 
moelectric  generator  Qgen  and  the  temperature  7}  seem  to  be 
independent  of  the  resistance  load,  the  two  contributions,  heat 
power  input  into  the  TEMs  and  heat  lost  through  the  assembly 
structures,  seem  to  be  also  independent  of  the  resistance  load. 

In  Tables  8  and  9  are  reported  the  different  components  of  the 
heat  flow  balance,  for  the  ‘open’  and  ‘closed’  circuit,  respectively. 
The  heat  losses  through  the  stainless  steel  screws  and  through 
the  thermal  insulation  between  the  TEMs  (filling  2/3  of  the  total 
area  of  the  generator)  are  both  significant  as  it  is  the  heat  ex¬ 
changed  to  the  ambient. 

The  ‘open’  and  ‘closed’  circuit  voltage  for  the  different  resistance 
loads,  grouped  for  different  values  of  the  heat  flow  into  the  TEMs 
(Tables  6  and  7)  are  shown  in  Fig.  7  as  a  function  of  the  tempera¬ 
ture.  The  ‘open’  circuit  voltage  increases  as  the  temperature  differ¬ 
ence,  AT,  increases,  and  is  well  interpolated  using  a  linear  function. 
The  same  occurs  for  the  ‘closed’  circuit  voltage  when  considering  a 
certain  resistance  load.  These  data  also  show  other  regular  behav¬ 
iour:  for  an  assigned  power  input  the  temperature  difference  in¬ 
creases  as  the  resistance  load  increases;  if  the  resistance  load 


Table  6 

Experimental  data  of  the  heat  power  input  for  the  different  data  set  (open  circuit). 


R2  =  22.39  n 


R:]  =  39.20  O 


&,  =  86.22  £2 


Qj  Qgen  QrEM  Qj  Qgen  QrEM  Qj  Qgen  QrEM  Qj  Qgen  QrEM 


107.9  101.8 

130.5  123.1 

167.9  157.3 

211.3  196.7 

237.3  219.7 


81.0  108.6 

97.8  131.0 

124.8  167.9 

155.9  212.0 

174.1  238.9 


102.6  81.6 

123.5  98.2 

157.4  124.9 

197.3  156.4 

221.3  175.3 


107.4  101.5 

129.4  122.1 

167.8  157.5 

211.2  196.9 

240.1  222.5 


80.7  107.6 

97.0  131.4 

125.1  167.1 

156.1  210.5 

176.2  241.3 


101.8  80 

124.1  98 

156.8  124 

196.3  155 

224.0  177 


Table  7 

Experimental  data  of  the  heat  power  input  for  the  different  data  set  (closed  circuit). 


R,  =  11.14n 


Qj  Qgen  QteM 

107.8  101.9  81.0 

130.7  123.4  98.1 

166.2  155.9  123.8 

212.0  197.7  156.8 

239.0  221.7  175.8 


R2  =  22.39  £2 

Qj  Qgen  QrEM 

109.5  103.6  82.4 

130.4  123.0  97.8 

168.2  157.8  125.3 

212.3  197.9  156.9 

240.8  223.4  177.1 


R3  =  39.20  O 
Q/Qgen 
108.0  102.3 

130.2  123.0 

165.8  155.7 

208.8  194.8 

242.9  225.4 


Qtem 

81.4 

97.8 

123.7 

154.4 

178.6 


R4  =  86.22  n 
Q(Qgen 

108.5  102.7 

130.8  123.6 

168.2  158.0 

208.7  194.8 

240.3  223.0 


Qtem 

81.7 

98.2 

125.4 

154.5 
176.7 
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Table  8 

Components  of  the  heat  flow  balance  for  the  circuit  closed  on  R3. 


Ql  (W)  QanvtW) 

108.0  5.7 

130.2  7.2 

165.8  10.1 

208.8  14.0 

242.9  17.5 


Qins  (W)  Qi,  (W) 

7.7  13.2 

9.3  15.8 

12.0  20.0 

153  25.0 

17.8  28.9 


Qtem  (W) 
81.4 
97.8 
123.7 
154.4 
178.6 


Table  9 

Components  of  the  heat  flow  balance  for  the  open  circuit. 


Qj  (W)  Qenv  (W)  Qins  (W)  Qj,  (W)  Qtem  (W) 


107.4  5.9 

129.4  7.3 

167.8  10.3 

211.2  14.3 

240.1  17.6 


7.7  13.1  80.7 

9.3  15.7  97.0 

12.2  20.3  125.1 

15.5  25.3  156.1 

17.7  28.6  176.2 


AT  (°C) 


Fig.  7.  Comparison  between  experimental  and  predicted  data. 


AT 


Fig.  8.  Reduced  temperature  difference  AT  versus  AT  and  RL. 


becomes  high,  the  ‘closed’  circuit  voltage  tends  to  the  correspond¬ 
ing  ‘open’  circuit  voltage. 

As  previously  emphasized,  for  each  value  of  the  heat  flow  into 
the  TEMs,  Qtem.  the  temperature  difference,  AT,  between  the  bot¬ 
tom  and  the  top  of  the  modules  increases  as  the  resistance  load  in¬ 
creases  (Fig.  7).  This  behaviour  can  be  understood  with  the  help  of 
the  theoretical  model,  by  comparing  the  different  contributions  in 
Eq.  (5).  The  Joule  contribution  is  always  smaller  than  those  due  to 


Peltier  and  Fourier  and  can  be  neglected.  For  the  same  heat  power 
input  into  the  modules,  if  the  load  resistance  increases  the  circulat¬ 
ing  current  decreases  and  thence  the  Peltier  contribution  de¬ 
creases,  whereas  the  Fourier  contribution  increases;  therefore  a 
higher  temperature  difference  is  experienced  (see  Appendix  A  for 
details). 

Fig.  7  shows  also  the  comparison  between  experimental  and 
predicted  data  for  the  ‘open’  and  ‘closed’  circuit  voltage  for  differ¬ 
ent  resistance  loads.  The  comparison  shows  that  the  results  given 
by  the  theoretical  model  are  in  very  good  agreement  with  the 
experimental  data.  The  temperature  difference,  AT,  is  measured, 
on  the  bottom,  on  the  assembly  surfaces  in  the  immediate  proxim¬ 
ity  of  the  TEMs,  and  on  the  top,  on  the  square  plate  of  aluminium 
over  the  TEMs. 

The  temperature  difference,  AT,  between  the  hot  and  cold  junc¬ 
tions  of  the  thermoelements  is  calculated  from  the  temperature 
difference,  AT,  with  the  procedure  explained  in  Section  3.  Fig.  8 
shows  the  trend  of  AT  versus  AT.  In  the  temperature  range  of 
the  experiments  a  linear  dependency  is  evident  of  AT  on  AT  for 
a  fixed  RL;  thus  the  predicted  ‘open’  and  ‘closed’  circuit  voltage 
are  linear  functions  of  the  temperature  difference  on  the  external 
surfaces  of  the  modules,  in  agreement  with  the  experimental  evi¬ 
dence  (see  Appendix  B  for  details). 

Fig.  9  shows  the  theoretical  and  experimental  power  output. 
The  data  distribution  is  typical  of  a  second  order  polynomial  of 
the  temperature  difference,  AT.  The  power  output,  P,  obtained  with 
the  theoretical  model  is  a  second  order  function  of  the  temperature 
difference,  AT  (Eq.  (4));  as  previously  observed  there  is  an  almost 
linear  correlation  between  AT  and  AT  (Fig.  8);  consequently,  the 
predicted  power  output  is  a  parabolic  function  of  AT,  in  agreement 
with  the  experimental  data  (see  Appendix  B). 

In  Fig.  9,  for  the  couples  of  resistance  R2  and  R3,  R,  and  R4  the 
power  output  shows  quite  overlapping  trends;  furthermore,  the 
measured  power  output  (Fig.  9)  is  greater  for  the  resistance  loads 
R2  and  R3.  This  behaviour  can  be  explained  by  Fig.  10,  where  the 
power  output  is  plotted  versus  the  resistance  load.  The  different 
curves  shown  in  Fig.  10  correspond  to  assigned  values  of  thermal 
power  input  into  the  TEMs,  Qtem.  This  is  obtained  as  the  average 
among  the  four  experimental  values  pertinent  to  the  five  power  in¬ 
puts  to  the  TEMs  (Tables  7  and  8).  The  different  curves  shown  in 
Fig.  10  have  a  maximum  between  R2  and  R3.  For  the  particular 
choice  of  the  load  resistances  R2  and  R3,  the  corresponding  values 
of  power  output  are  quite  equal.  For  the  couple  Ri  and  R4,  a  differ¬ 
ence  occurs.  This  difference  increases  when  QTEm  increases  and  can 
attributed  to  the  different  AT  evident  in  Fig.  7  for  a  same  Qtem.  as 
discussed  in  Appendix  A.  The  maximum  power  output  is  achieved 
for  a  load  resistance  of  about  30  Cl,  whereas  the  inner  resistance  is 
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27  a  This  is  in  accordance  with  Freunek  et  al.  [23]  and  implies  that 
the  load  resistance  enabling  the  generator  to  operate  at  the  maxi¬ 
mum  output  power  is  equal  to  the  internal  one  added  of  a  term 
function  of  AT/ AT. 

The  comparison  of  predictions  and  measurements  based  on 
electric  power  data  (Figs.  7-10)  is  particularly  good.  When  moving 
to  the  heat  exchange,  the  accordance  is  no  longer  so  good.  In  Fig.  1 1 
the  calculated  heat  power  into  the  modules,  Qtem,  and  into  the 
generator,  Qgen,  are  compared  with  the  corresponding  experimen¬ 
tal  data.  The  agreement  is  good  in  terms  of  tendency,  even  if  quan¬ 
titatively  poor.  The  calculated  values  are  always  lower  than  those 
measured.  This  difference  is  certainly  due  to  a  basic  difficulty  in 
assessing  the  thermal  resistances  in  the  experimental  setup.  Actu¬ 
ally,  these  are  calculated  referring  to  the  available  schemes  for  the 
cooling  of  electronic  systems  [24]  and  seem  to  be  lower  than 
expected. 

In  Fig.  12  are  shown  the  experimental  values  of  the  net  effi¬ 
ciency  of  the  TEMs,  £tem,  as  a  function  of  the  temperature.  For  each 
resistance  load,  the  data  distribution  show  a  linear  trend.  In  Fig.  12 
are  also  shown  the  values  of  the  net  efficiency  calculated  by  means 
of  the  theoretical  model.  The  predicted  values  are  higher  than  the 
experimental  data  (Fig.  12),  this  being  in  accordance  with  that  pre¬ 
viously  observed.  Referring  to  Eq.  (16),  for  the  power  output  (the 
numerator  in  Eq.  (16))  in  Fig.  9  there  is  a  good  agreement  between 
theoretical  and  experimental  data;  for  the  heat  flow  (the  denomi¬ 
nator  in  Eq.  (16))  in  Fig.  11  the  predictions  are  lower  than  the  mea- 


Qgen,exp(W),  QTEM,exp(W) 


Table  10 

Experimental  efficiency  (%). 

Rl  =  11.14  Q  RL  =  22.39  fi  Rt  =  39.20  Cl  RL  =  86.22  Cl 

(fi)net  (c)ov  (£)net  (£)ov  (£)net  (£)ov  (£)net  (£)ov 

28 
34 
43 
52 
59 


surements.  This  makes  the  predicted  efficiency  higher  than  the 
measured  one. 

In  Table  10  are  reported  the  net  and  the  overall  efficiency  of  the 
thermoelectric  generator.  The  net  efficiency  does  not  exceed  1% 
while  the  overall  efficiency  reaches  0.7%.  The  lower  overall  effi¬ 
ciency  is  basically  due  to  the  presence  of  the  thermal  insulation 
among  the  modules.  Furthermore,  the  measurement  of  the  tem¬ 
perature  on  the  hot  side  of  the  modules  requires  a  thick  layer  of 
thermal  grease  between  the  top  plate  and  the  TEMs.  This  thick 
layer  penalizes  both  the  net  and  the  overall  efficiencies.  Reducing 
this  thickness  in  the  model,  the  net  efficiency  increases  to  about 
20%  and  the  overall  efficiency  to  about  25%. 

In  this  work  thermoelectric  cooling  modules  (TECs)  are  used  as 
thermoelectric  generators  (TECs).  In  order  to  discuss  this  choice, 
TECs  and  TECs  are  compared  in  terms  of  power  output  and  conver¬ 
sion  efficiency  for  the  same  condition  (TH  =  90  °C  and  Tc  =  30  °C). 
The  thermoelectric  modules  are  reversible  and  can  be  used  both 
as  TECs  and  as  TEGs.  However,  different  applications  require  spe¬ 
cific  arrangements. 

In  terms  of  themo-mechanical  features,  since  the  TEGs  are  sub¬ 
mitted  to  higher  temperature  differences,  the  electrical  wires  are 
soldered  on  the  cold  side  of  the  module  and  the  insulating  materi¬ 
als  of  the  wires  are  rated  for  250  °C.  Since  the  TECs  are  submitted 
to  higher  circulating  currents,  the  electrical  wires  of  the  TECs  are 
thicker  as  their  insulating  materials,  rated  for  90  °C.  Finally,  the 
price  of  the  TECs  is  sensibly  lower  than  that  of  the  TEGs  [3], 

In  terms  of  thermo-electrical  features,  in  Fig.  13  is  shown  a 
comparison  of  different  commercial  TECs  and  TEGs  (Table  1 1 )  in 
terms  of  power  output  and  overall  conversion  efficiency.  The  TEC 
used  in  this  experiment  gives  performances  sensibly  lower  than 
the  other  modules.  The  maximum  power  produced  by  the  remain¬ 
ing  modules  is  not  so  far  and  it  is  evident  that  in  the  present  test 
conditions,  if  the  quality  of  the  semiconductors  is  high  enough,  a 
TEC  can  produce  as  a  TEG.  It  can  be  concluded  that  TECs  are  an 
interesting  alternative  for  power  generation  in  the  “low-heat” 
region  [3],  Obviously,  due  to  the  thermo-mechanical  features,  if 


Fig.  11.  Experimental  and  predicted  power  heat  ir 
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Fig.  13.  Comparison  of  power  output  (a)  and  overall  efficiency  (b)  for  different  TECs 
and  TEGs  in  the  same  conditions  (AT  =  60  °C,  Tc=  30  °C). 


This  experimental  work  gives  useful  information  to  increase  the 
effectiveness  of  devices  for  thermoelectric  conversion. 

It  may  be  generally  concluded  that  thermoelectric  conversion  is 
a  promising  way  to  the  recovery  of  waste  heat,  because  there  is  sig¬ 
nificant  room  for  improvement:  the  use  of  suitable  modules  for  the 
Seebeck  mode,  a  greater  care  of  the  thermal  contact  with  the  hot 
and  the  cold  source  and  a  greater  density  of  the  modules  in  the 
generator.  Furthermore,  the  conversion  efficiency  reached  in  the 
range  of  the  temperature  difference  of  the  experiment  is  rather 
low,  so  the  trend  of  the  gathered  data  allows  the  behaviour  of  a 
generator  with  higher  heat  power  input  to  be  predicted. 
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Appendix  A 

Eq.  (5)  can  be  simplified  neglecting  the  thermal  resistance  of  the 
ceramic  plates  and  the  Joule  contribution,  always  very  small.  Then 
Eq.  (5)  can  be  rewritten  as: 

QtEm  =  mN  [X'  (T'H  -  Tc)  +  (a,  -  <Xn) 7V]  (A.1) 

Eq.  (A.1 )  can  be  rewritten  to  put  in  evidence  the  contribution  due  to 

VH: 

^  +  K'T'C  JT  +  I(ap  -  a„)]  T„  (A.2) 

Since  the  bottom  temperature  Tc  is  quite  constant  (Tables  5  and  6) 
and  the  thermal  conductance  Kc  is  very  high,  the  temperature  T'c  is 
also  quite  constant,  as  is  Qti:m  (Tables  7  and  8).  Thence,  the  left  term 
of  Eq.  (A.2)  is  constant  and  when  the  circulating  current  I  increases 
the  temperature  TH  decreases. 


Table  11 

Thermocouple  properties  for  TEC  and  TEG  module. 


Module  N  Size  (mm)  a  K  R  Z 

(pVK-1)  (mWr1)  (mil)  (K  1 ) 


Test  module  127  40x40x3.8  343  3.53 

TEC  127  couples  127  40x40x4.0  433  3.89 

TEG  254  couples  254  40x40x3.8  441  2.92 

TEG  127  couples  127  40  x  40  x  3.6  380  3.42 


23.7  0.00141 

21.8  0.00210 

31.2  0.00213 

18.4  0.00229 


higher  temperature  differences  are  available,  TEGs  become  the 
right  choice. 


6.  Concluding  remarks 


Appendix  B 

In  Fig.  8  is  shown  the  temperature  difference  on  the  junctions, 
AT”,  as  a  function  of  that  outside  the  ceramic  plates,  AT,  and  of 
the  resistance  load,  RL.  The  dependency  is  linear  and  in  first 
approximation,  can  be  expressed  as: 

AT'  =  C(Rl)AT  (B.l) 

where  the  constant  C  is  a  function  of  Rl.  Consequently,  Eqs.  (2),  (4), 
and  (8)  can  be  expressed  as  a  function  of  AT.  For  example  the 
‘closed’  circuit  voltage  becomes: 

Vcc  =  (RTEMRL+Rt)  mN^P  -  «n)C(RL)AT  (B.2) 


An  experimental  investigation  has  been  carried  out  of  the  per¬ 
formance  of  multiple  conventional  Peltier  devices  used  in  the  See¬ 
beck  mode  as  a  thermoelectric  generator.  This  was  tested  within  a 
purpose-built  heat  transfer  closed-loop. 

The  good  agreement  between  the  experimental  and  calculated 
values  of  the  different  parameters  allows  the  use  of  the  theoretical 
model  as  a  tool  of  analysis  of  the  generator  performance.  The  low 
net  efficiency  is  attributed  to  the  heat  transfer  not  being  optimal 
and  to  the  low  quality  of  the  modules.  The  commercial  modules  used 
in  the  experiment  have  a  low  figure  of  merit  (Z  =1.41  x  10  3  K  '), 
while  for  an  efficient  power  generation  higher  values  are  preferred. 

The  comparison  between  the  net  and  the  overall  efficiency 
highlights  the  negative  role  of  the  insulating  layer  and  allows  the 
reduction  of  the  generator  efficiency  to  be  estimated. 
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